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Within the framework of Dyson-Schwinger equations (DSEs), we discuss the equation of state 
(EOS) and quark number densities of 2+1 flavors, that is to say, u, d, and s quarks. The chemical 
equilibrium and electric charge neutrality conditions are used to constrain the chemical potential of 
different quarks. The EOS in the cases of 2 flavors and 2+1 flavors are discussed, and the quark 
number densities, the pressure, and energy density per baryon are also studied. The results show 
that there is a critical chemical potential for each flavor of quark, at which the quark number density 
turns to nonzero from 0; and furthermore, the system with 2+1 flavors of quarks is more stable than 
that with 2 flavors in the system. These discussion may provide some useful information to some 
research fields, such as the studies related to the QCD phase transitions or compact stars. 
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I. INTRODUCTION 

Quantum Chromodynamics (QCD) is commonly ac¬ 
cepted as the basic theory of strong interaction. Because 
of the asymptotic freedom character of QCD, the high 
energy processes could be pictured hadrons as weakly in¬ 
teracting quarks and gluons. However, the interaction 
is so strong that the perturbation theory is invalid at 
energies below a few hundred MeV, and the picture is 
inadequate. Further more, confinement and dynamical 
chiral symmetry breaking (DCSB) are two fundamental 
properties of QCD in the region of low energy. On the 
one hand, the color confinement is not understand very 
well until now, it remains only an empirical fact that col¬ 
ored object have not been observed. On the other hand, 
DCSB is better understood, which explains the origin 
of constituent-quark masses and underlies the success of 
chiral effective field theory. 

The study of QCD at nonzero temperature T and fi¬ 
nite baryon density pB is one of the most interesting 
topics in contemporary physics studies. It is believed 
that a phase transition happens in strongly interaction 
matter from a confined hadronic system at low temper¬ 
atures and low densities to a deconfined quark gluon 
plasma (QGP) phase at high temperatures and/or den¬ 
sities. There are several famous laboratories and ex¬ 
perimental facilities which are related with thermal and 
dense QCD, such as the Relativistic Heavy Ion Collider 
(RHIC), the Large Hadron Collider (LHC), and the Fa¬ 
cility for Antiproton and Ion Research (FAIR). RHIC 
and LHC are mainly concentrated on high temperature 
and low density physics, while FAIR focus on low tem¬ 
perature and high density plasma physics. It is known 
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that the Equation of state (EOS) of QCD for cold and 
nonvanishing quark chemical potential p plays a key role 
in the study of compact stars [1-4]. Since the perturba¬ 
tion QCD is invalid in the low-energy region, it is very 
difficult to obtain a reliable EOS from the first princi¬ 
ples of QCD. A lot of models of neutron stars or quark 
stars have been studied [5-7]. Many important studies 
have been done in Refs. [ 8 , 9], and as data of masses 
and radii of neutron stars have accumulated, the EOS 
that best fits the observation data has been strictly con¬ 
strained. The calculations of the EOS [10-15] and the 
transition temperature [16-18] based on the Lattice QCD 
can be performed with merely realistic quark mass spec¬ 
trum. However, Lattice QCD confront sign problem, 
which refers to the difficulty of evaluating the integral 
of a highly oscillatory function numerically, when deal¬ 
ing with finite chemical potential. Consequently, nowa¬ 
days various effective models are often and in some sense 
inevitably used to study related issues phenomenologi¬ 
cally, such as the Nambu-Jona-Lasinio (NJL) model [19- 
31], Lattice QCD [32-34], the Dyson-Schwinger equa¬ 
tions (DSEs) [35-44], and the Quantum Electrodynam¬ 
ics in (2+1) dimensions (QED 3 ) [35, 45-48]. This study 
will mostly focus on the EOS at zero temperature and 
nonzero chemical potential, especially the region of high 
chemical potential, within the framework of DSEs. One 
of our main goal is to compare the stability of the dense 
system between 2 flavors (only the two lightest quarks 
are included, namely, u and d quarks) and 2+1 flavors 
(not only u and d quarks, but the slightly heavier s quark 
is also included) cases. 

In this paper, the quark number density in the case 
of 2+1 flavors is calculated by constrains of the chemi¬ 
cal equilibrium and electric charge neutrality conditions 
within the framework of DSEs. The Qin-Chang effective 
gluon propagator model [49], which has been proved to 
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be successive in hadron physics and thermal QCD recent 
years, is adopted to study the EOS in both 2 flavors and 
2+1 flavors cases. The following of this paper is organized 
in such a way: In Sec. II, we give a brief introduction to 
the DSEs at zero temperature as well as zero chemical 
potential first, and then fit the quark propagator using 
the formula with 3 complex conjugate poles. In Sec. Ill 
we discuss in detail the quark number densities at T = 0 
algebraically, and show numerical results of the pressure 
in cases of 2 flavors and 2+1 flavors. Finally, we give a 
brief discussion and summary in Sec. IV. 


II. DYSON SCHWINGER EQUATIONS AT 
ZERO TEMPERATURE AND ZERO CHEMICAL 
POTENTIAL 


In this section, we briefly introduce the Dyson 
Schwinger equations (DSEs) formalism, which are cou¬ 
pled integral equations relating the Green’s functions for 
the theory to each other [35-38, 40]. The zero temper¬ 
ature and zero quark chemical potential version of the 
quark propagator DSE, which is the most important one 
among all the DSEs, reads 

S{p)~^ = Z2{if + Z^m) 

^g^ZiF J Y'^tiS{q)Tl{p,q)Df,^{p - q), ( 1 ) 

where S{p) is the dressed quark propagator, Z 2 is the 
field-strength renormalization constant, Z^ is the mass 
renormalization constant where m is the current quark 
mass, g is the coupling constant, Zip is the quark-gluon- 
vertex renormalization constant, := J is a sym¬ 
bol that represents a Poincare invariant regularization of 
the four-dimensional Euclidean integral, A“ is the Gell- 
Mann matrices, r“(p, g) is the quark gluon vertex, and 
D^v{p — q) is the dressed gluon propagator. The general 
form of the quark propagator at zero temperature and 
zero chemical potential reads 

S{p)-^ =iM{P^)+B{p^), (2) 

where A{p^) and B{p^) are scalar functions of p^. The 
renormalization condition is 

AiC) = I, (3) 

B{0 = m, (4) 

at sufficiently large spacelike [50]. We choose the 
renormalization point to be C = 19 GeV in this work. 
Eq. (1) is the exact result from first principles of QCD, 
but we can not solve it directly unless concrete trunca¬ 
tions are performed. Rainbow truncation is used in this 
work, which means a bare vertex is adopted 

r“(P>9) = (5) 


and the Qin-Chang gluon propagator model [49] is spec¬ 
ified by a choice for the effective interaction in Landau 
gauge. 


with 
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the transverse projection operator, and 
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(7) 
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(9) 

( 10 ) 


where the parameters are jm = 12/(33 — 2Nf) with 
Nf = A, T = - 1 with e = 2.71828, = 0.234 

GeV, and mt = 0.5 GeV, which determined by behavior 
of ultraviolet of the gluon propagator. The remaining 
parameters D and lo are fitted by the mass and decay 
constant of pion, we choose w = 0.7, Duj = (0.8GeV)^ 
in this work as in Ref. [49]. Moreover, in our study the 
isospin symmetry is preserved in the strong interaction 
section, and the current quark mass m„ = = 3.4 

MeV, nis = 82 MeV is used. The isospin symmetry will 
be broken when the electroweak interactions is included, 
which will be introduced in Sec. Ill of this paper. 

Now substitute Eqs. (2), (5), and ( 6 ) into Eq. (1), and 
project Eq.(l) on the two tensor structures Ti = ip and 
T 2 = 1, one can obtain the following equations for the 
two dressing functions A{p^) and B{p^), 


A{p^) — Z 2 -\- 


Qik^) 


A{q^) 


3p2 


k2 q2 {q^) + B'^ [q^) 


X [p- q + 2 


{k-p){k ■ q) 


fc 2 


B{p'^) = Z 4 m + 4 / 

J a 


G{k^) 


B{q^) 


k2 {q^) + B"^ [q^) 


( 11 ) 
■ ( 12 ) 


where k^^ = p^ — g^.The coupled Eqs. (11) and (12) can 
then be solved by direct iteration. The DSE solution 
for the dressed quark propagator can be well fitted with 
3 pairs of complex-conjugate poles with the representa¬ 
tion [51-54] 


•5(P) = V 


n—1 


Zn 

ip + TO„ 



(13) 


Recently, this formula is generalized to study parton dis¬ 
tribution amplitude of mesons [55, 56]. We fit u, d and 
s quarks respectively, where the following requirement, 
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that the quark propagator in the region of ultraviolet 
should tend to the free quark propagator, is employed 

^(zfc+4) = l. (14) 

fe=i 

Corresponding parameters for u, d quarks are 

mi = 349 + llli MeV, = 0.3874 + 0.618i, 
m 2 =-1224 + 272i MeV, Za = 0.0912 + 0.176i, (15) 
m 3 = 1356 + &Mi MeV, Z 3 = 0.0214 - 0.097h 

and for s quark, 

mi = 714 + 215* MeV, zi = 0.4060 + 0.689*, 

m 2 =—1863 — 96* MeV, 2:2 = 0.0987 + 0.0952*, (16) 

m 3 = 1429 + 893* MeV, Z 3 = -0.0047 - 0.093*. 



FIG. 1. Quark number density as a function of for u, d 
quarks, where T is taken to be 0. 


III. QUARK NUMBER DENSITY AND EOS IN 

Nf ^2 AND Nf = 2 + 1 CASES 

The solution of the quark propagator at T = 0, ^ = 0 
case can be generalized to T 7 ^ 0 , /* 7 ^ 0 case by the 
following replacement [57] 

p'^ ^ Cdn = LUn + (17) 

which is widely used within rainbow truncation of DSEs 
in thermal and dense QCD, and a;„ = (2n + l) 7 rT is 
the Matsubara frequencies. It has been proved [57] that 
based on the rainbow truncation of the DSEs and two 
assumptions: (1) The full inverse quark propagator at 
finite chemical potential is analytic in the neighborhood 
of p= 0; (2) The effective gluon propagator is /* indepen¬ 
dent, the quark propagator at nonzero T and p can be 
obtained by replacement —>■ Wn = + *f*- Hence, the 
quark propagator at nonzero T and /* is 



FIG. 2. Quark number density as a function of p for s quark, 
where T is taken to be 0. 

The summation of infinite Matsubara frequencies can be 
evaluated using the method of contour integral [60] 


S{p,COn) = XI ( “ 


Zk 


^ V *7 • P + + ruk 


i+ ■p + i^iUJn + ml 


(18) 


The well-known formula for the quark number density 
is [58, 59] 


+ 00 


p{p,T) =-N,Nj T 


d?p 

(^ 


tr[S'(p,w„) 74 ]. (19) 


Substituting Eq.(18) into Eq.(19), and then performing 
the trace, one obtains 
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= AiN^Nf 




k^l 


{2tt) 


Zk{uJn + ip) 


J? + {ujn + ipY + ‘ 


+ *p) 


p2 + {ujn + ipY + ml^ 


(20) 


+ 00 

T X ■^(^0 = iw„) 

n— — oc 


= -9 X 


Re(p(j)>0 


/(po) tanh ^ + /(-po) tanh ^ 


( 21 ) 


where po is regarded as a complex variable and Res[- • •] 
means the residue of the function. Taking the limit T -+ 
0 of p{p,T), one obtains 


p{p,T = 0) 

= “ V ~ 

( 22 ) 

where = |R.e(mfc)| and Ck, dk are defined by m^ = 

Ck + dff^h. 

The quark number densities dependence on p at T = 0 
for u, d and s quarks are displayed in Fig. 1 and Fig. 2 


























4 


respectively, u, d quarks have the same quark number 
density dependence on their own chemical potential since 
the isospin symmetry is used and then their current quark 
mass are the same. We can see from Fig. 1 that there is 
a critical chemical potential = 349 MeV, at which the 
quark number density of u, d quarks turn to nonzero from 
exact zero, while Fig. 2 shows that the quark number den¬ 
sity of s quark turns to nonzero at Hc = 714 MeV.^ This 
behavior agrees qualitatively with the general conclusion 
of [61]. 

According to the electroweak theory of the Standard 
Model, 3 flavors of light quarks can transform with each 
other by the following reactions [62]. 


d ^ u + e~ + Vg, 

(23) 

u + e~ d + Vg, 

(24) 

s ^ u + e~ + Vg, 

(25) 

u + e~ —>■ s + Z/e, 

(26) 

u + d u + s . 

(27) 


Since different quarks are coupled with each other, we 
can then employ the chemical equilibrium and electric 
charge neutrality conditions to constrain different chem¬ 
ical potentials, that is to say, /x„, jig and The 
conditions read 


fj-d = fj-u + fj^e, (28) 

/ig — 4“ Me; (29) 

fpti - ^Pd - \ps - Me = 0. (30) 

There is only one independent chemical potential due to 
constraint equations Eqs. (28), (29) and (30), we choose 
/r„ in this work. As discussed above, the quark number 
densities of u, d and s quarks at T = 0 and finite p are 
showed in Fig. 1 and Fig. 2, and the particle number 
density of electron at T = 0 is an explicit function of Me; 
namely [60] 


PeiPe) = (31) 

For definite chemical potential of quark, the EOS of QCD 
at T = 0 reads [59, 63] 

P{p) = P{p = 0)+ r dp'pip'). (32) 

Jo 



FIG. 3. The particle number density dependence of u, d, 
and s quarks, as well as electron, on the chemical potential 
of u quark, by the constrain of the chemical equilibrium and 
electric charge neutrality conditions. 



FIG. 4. The pressure dependence of u, d, and s quarks, as 
well as electron, on the chemical potential of u quark, by 
the constrain of the chemical equilibrium and electric charge 
neutrality conditions. 


And in the case of 2+1 flavors, where s quark is included, 
the pressure is 

P{pu) = Pu{pu) + Pd{pd) + Ps{ps) + Pe{pe)^ ( 34 ) 

where pd, Ps, Pe are all related to pu- The baryon num¬ 
ber density in this case is then defined by 

Of I . _ Pu{Pu) + Pd{Pd) + Ps{Ps) 

Pb {Pu) = - 5 -■ ( 35 ) 


In the case of 2 flavors, the EOS is just take Nf = 2 for 
the two light quarks, u and d, and the baryon number 
density is defined by 

2/ _ Puipu) + Pd{pd) /'qq'i 

Pb = -5-■ (33) 


^ In this work, the value of fic depends on the real part of the first 
mass pole. 


Now we plot the particle number density dependence 
and the pressure dependence of u, d, and s quarks, as 
well as electron, on u quark chemical potential pu in 
Fig. 3 and Fig. 4 respectively. The pressure depen¬ 
dence of the system on baryon number density in cases 
of 2 flavors and 2+1 flavors are shown in Fig. 5. The 
quark number density of u and d quarks turn to nonzero 
at Pu = 349 MeV with pu increasing. In the area of 
Pu ^ 610 MeV, the quark number density of s quark is 
exact 0. It indicates that at this range the EOS and en¬ 
ergy density of dense system has no difference between 
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PflCMeV^) 



P„ (MeV) 


FIG. 5. The pressure dependence on baryon number density 
for Nf = 2 and Nj = 3 cases. 



Pb 

FIG. 6. The dependence of energy density per baryon on 
baryon number density for Nf — 2 and Nf = 3 cases. 


2 flavors and 3 flavors cases. Combining with Fig. 4, 
we can see that the pressure of u, d, and s quarks are 
Pu = 9.74 X 10® MeV'‘, Pd = 2.65 x 10^ MeV'‘, and 
Ps = 0 respectively. We deduced from Fig. 3 and Fig. 
4 that the s quark has no effect on the pressure of the 
system when ;< 610 MeV, and the critical chemical 
potential of u quark, /r„ = 610 MeV, corresponds to 
baryon number density ps = 1.027 x 10^ MeV® and the 
pressure of system Pt = 3.62 x 10® MeV^. Such conclu¬ 
sion is just depicted in Fig. 5 that the two curves are 
indistinguishable when ps ^ 1.027 x 10^ MeV®. Two 
curves begin to split at ps = 1.027 x 10^ MeV®, namely 
p-u = 610 MeV, and the pressure in the case of 2-1-1 fla¬ 
vors is smaller than that in the case of 2 flavors at the 
area of pb ^ 1.027 x 10^ MeV®. 

In the limit of large pu, or in other words, large pm Pd, 
Ps simultaneously, the behavior of pressure should tend 
to NcPq /like the free quark gas, where pq means 
the chemical potential for it, d, or s quarks. This can 
be easily understand since the quarks have the feature 
of asymptotic freedom. The pressure compared with free 
quark gas in 2-1-1 flavors and 2 flavors is depicted in Fig. 7. 
We can see from Fig. 7 that the pressure of 2-|-l flavors 


FIG. 7. The pressure compared with free quark gas for Nf — 2 
and Nf = 3 cases. 

case tends to that of the free quark gas more slowly as the 
chemical potential increases. Compared with the work by 
Fraga et al, our results imply that the interaction in the 
case of 2-1-1 flavors is stronger than 2 flavors case and that 
in the perturbative EOS [64, 65]. It is consistent with the 
view of strong coupled quark-gluon plasma, even in the 
case of large chemical potential. 

The relation between the energy density and the pres¬ 
sure of the corresponding system is [66, 67] 

£ = -PP^PiP,, (36) 

i 

where pi and pi represents the chemical potential and 
the particle number density for each component in the 
system. We plot the dependence of energy density per 
baryon with the baryon number density in Fig. 6. It is 
shown that the energy density per baryon is much smaller 
for the 2-1-1 flavors than that for the 2 flavors case when 
the s quark number density is nonzero. Our calculations 
imply that the system with 2-1-1 flavors of quarks is more 
stable than the case with 2 flavors of quarks. 

IV. DISCUSSION AND SUMMARY 

In this paper, the equation of state of 2-1-1 flavors of 
quarks, that is to say, it, d, and s, is studied within the 
framework of Dyson-Schwinger equations. We solve the 
DSE of 2-1-1 flavors quark propagator under the famous 
Rainbow-Ladder truncation, which preserves axial vec¬ 
tor current conservation. An effective gluon propagator 
model, namely the Qin-Chang model [49], is adopted. It 
has been proved to describe hadron physics very well, 
especially for ground states. The gluon model generates 
DCSB, which is an important feature of low energy QCD. 
And it is infrared finite, which is consistent with lattice 
QCD’s results [49, 68, 69]. The parameters of this model, 
i.e., iriu, md, D, and w are determined by fitting the 
masses and decay constants of pion and kaon. It is worth 
noting that these parameters are fixed during our study. 
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The three complex-conjugate poles formula is used to fit 
the quark propagator at zero temperature and zero quark 
chemical potential. The quark propagator at finite tem¬ 
perature and finite chemical potential is generalized via 
the standard way in finite temperature field theory us¬ 
ing the Matsubara frequencies. The chemical potential 
dependence of quark number densities are calculated at 
T = 0. As Ref. [61] pointed out, the existence of some 
singularity at the point /x = /ic, T = 0 is a robust and 
model independent prediction, and the critical chemical 
potential for the lightest quarks (u and d) is estimated 
to be /X = /ic « TOjv — 16 MeV. In our study, the value of 
critical chemical potential is located at /x = 349 MeV, in 
which the quark number density turns to nonzero from 0 
with the u quark chemical potential increasing. It is qual¬ 
itatively consistent with the general conclusion in [61]. 
As the quark chemical potential increases, the pressure 
of our result tends more slowly to the one of the free 
quark gas. This suggests that the interaction is stronger 
than the corresponding one in perturbation EOS, which 
is in consistent with the view of strongly coupled QGP, 
even in the case of large chemical potential. 

It is well known that different quarks can transform 
with each other by electroweak interactions, the chemi¬ 
cal equilibrium and electric charge neutrality conditions 
are used to constrain the different chemical potentials, 
namely /x„, Ms and /Xe- We choose Mm as the inde¬ 
pendent variable, in other words, Md, Ms and Me are all 
functions of Mu- In the cases of 2 flavors and 2-1-1 flavors, 
the quark number densities and pressure dependence on 


the chemical potential of u quark are calculated, respec¬ 
tively. It is displayed in Fig. 3 that ps is exactly zero 
when /Xu < 610 MeV, so that the case with 2-1-1 flavors 
has no difference with the case of 2 flavors as showed in 
Fig. 5. While in the region of pB > 1.027 x 10^ MeV^, 
the pressure in the case of 2-1-1 flavors is larger than the 
case with 2 flavors. Subsequently, the dependence of en¬ 
ergy density per baryon on baryon number density for 
Nf = 2 and Nf = 3 cases is studied. Fig. 6 shows that 
the energy density per baryon with 2-1-1 flavors of quarks 
is smaller than that with 2 flavors. It implies that system 
with 2-1-1 flavors of quarks is more stable in the region 
where ps > 1.027 x 10^ MeV^. Our results imply that 
the interaction in the case of three flavor is stronger than 
two flavors case and perturbation EOS. 

Finally, we stress the two key points in this paper. 
Firstly, we generalize the EOSs in SUf{2) to that of 
SUf{3). Furthermore, We firstly combine DSEs with 
chemical equilibrium and electric charge neutrality con¬ 
ditions to study EOS in dense QCD. These studies may 
provide some useful information to some research fields, 
such as the studies related to the QCD phase transitions 
or compact stars. 
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